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Abstract 

In the present article, we study the numerical approximation of a system of Hamilton- Jacobi 
and transport equations arising in geometrical optics. We consider a semi-Lagrangian scheme. We 
prove the well posedness of the discrete problem and the convergence of the approximated solution 
toward the viscosity-measure valued solution of the exact problem. 
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1 Introduction 



In this article, we consider the following system 

dtu + H{x,t,Vu) = 0, 
dtm -I- V • (a(x, Vu) m) = , 
u{x,Q) = uo{x), ?n(0) = mo , 



in M'^ X (0,r], 
in M'^ X (0,r], 
in M'^ , 



(1.1) 
(1.2) 
(1.3) 



of an Hamilton-Jacobi type equation (HJ) and a continuity equation (CE) describing the transport of 
the conserved measure mg. Even if the vector field a can be smooth (in the simplest and of reference 
case, a{x,p) = p), the scalar function n as a solution of (1.1) is intended in the viscosity sense. 
Therefore, Vu is at most L°°, even for smooth initial data uq, and the regularity that we can expect 
for the vector field a(-, Vu) is the very weak regularity 



x[0,T])r 



(1.4) 



As a consequence, despite the fact that the (HJ) equation can be solved independently from (CE), 
system (1.1)-(1.2) leads to some interesting mathematical issues as well as numerical challenges. 

More specifically, in order to obtain global solutions, a well adapted notion of solution for (1.2) 
under hypothesis (1.4) is the notion of measure solutions introduced by Poupaud-Rascle [16]. Indeed, 



the latter perfectly makes sense if the vector field a(-,Vn) satisfies a one-sided Lipschitz condition 
(OSL), (see (2.3) below), and this condition can be obtained by the semiconcavity of the viscosity 
solution of (1.1), at least in the reference case a{x,p) = p, or in the one dimensional case for a class 
of functions a. 

Existence and uniqueness results for problem (1.1)-(1.2)-(1.3) in the framework of viscosity-measure 
solutions have been given for example in [2, 19] and in the different framework of the one dimensional 
viscosity-duality solutions in [12]. Therefore, our goal here is not to refine these previous results but 
to construct consistent numerical approximations. 

The considered scheme is based on a semi-lagrangian discretization of (1.1)-(1.2) for the time ap- 
proximation coupled with a finite element discretization for the space variable. It requires a convex 
hamiltonian H. Since the semiconcavity of the initial data uq is conserved by the scheme, it is sufficient 
to require only a weak (OSL) condition at the discrete level (automatically satisfied in the reference 
case a{x,p) = p) without semiconcavity requirement for the viscosity solution u. 

The stability of the Filippov characteristics and of the corresponding measure solution of the trans- 
port equation (CE) will be the key tool to prove the convergence of the numerical schemes. As 
a by-product, we obtain of course a new existence and uniqueness result for the viscosity-measure 
solution of (1.1)-(1.2)-(1.3). 

Regarding the applications of our numerical approximations of (1.1)-(1.2)-(1.3), it is worth to recall 
that these types of systems arise for example in the semi-classical limit for the Schrodinger equation 
[4] and of the spinless Bethe-Salpeter equation (the relativistic Schrodinger equation, [3]), or in the 
high frequency approximation of the Helmholtz equation. In these cases, the hamiltonian H and the 
vector field a take the forms 

H{x,t,p) = -\p\^ + V{x,t) , a{x,p)=p, {Schrodinger and Helmholtz equation), (1.5) 



where V{x,t) is the potential (see [2, 12, 19] for a rapid derivation of (1.5) and (1.6)). 

Coupling between first or second order Hamilton-Jacobi equations and transport equations has been 
also recently considered in the framework of Mean Field Games theory [13] . In this case however, the 
system is given by a backward Hamilton-Jacobi equation and a forward transport equation with 
initial-terminal conditions and coupling terms in both equations. The resulting system turns out to 
be completely different from the one considered here. 

The paper is organized as follows. Section 2 is devoted to the preliminary definitions and known 
results concerning system (1.1)-(1.2)-(1.3). In Section 3 we construct a semi-lagrangian scheme for 
the time discretization of system (1.1)-(1.2)-(1.3) and we prove its convergence. The corresponding 
fully discrete scheme is given and analyzed in Section 4. Finally, the appendix is devoted to the proof 
of some technical lemmas for the reader convenience. 

Throughout the paper, we will denote by Cq (M'^) (resp. C^(M'^)) the space of continuous functions 




{Bethe-Salpeter equation), (1.6) 
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which tend to at mfinity (resp. with compact support); by a standard mollifier, i.e. /5e(x) = ^p(|), 
p € C^(M'^), p > and p{x)dx = 1; by * the convolution with respect the space variable and by 
C any numerical constant that can vary from line to line in the computations. 

2 Preliminaries : the viscosity-measure solutions 

As mentioned in the introduction, a solution of (1.1)-(1.2) is intended in the viscosity sense for (1.1), 
while in the sense of Poupaud-Rascle [16] for (1.2). Concerning the definition of viscosity solutions, 
we refer to the pioneering articles [6, 5]. Here, for the reader's convenience, we recall the usual 
assumptions on H and the consequent general existence and uniqueness result for (1.1) that we shall 
use in the sequel. 

Let us define Qt '■= M"^ x [0,T], T > 0, d > 1, and let the hamiltonian H satisfy the following: 
(Hi) H is uniformly continuous on Qt x B{0, R) for any R > ; 
(H2) H{x,t,0) is uniformly bounded : supg^ \H{x,t,0)\ = M < +00 ; 

(H3) there exists rj > s.t. : \H{x,t,p) - H{y,t,p)\ < 77(1 + \p\)\x - y\ , V t G [0,T] , V x,y,p G . 

Theorem 2.1 ([6]). Under hypothesis (Hi), (H2) and (H3), if in addition the initial data uq belongs 
to (W'^'°°nBUC){R'^), there exists a unique viscosity solution u G {W^^°° r\BUC){QT) of (1.1), (1.3). 

As stated in Theorem 2.1, the expected regularity for is L°°{Qt)- Therefore, the characteris- 
tics X{t;x) associated to the conservative transport equation (1.2) cannot be defined as classical or 
distributional solutions of the system below 

dtX(t; x) = a(X(t; x),S/u(X(t] x),t)) 

(2.1) 

X{0;x)=x, xeR'^, 

but have to be understood in a generalized sense. Once the flow X{t;x) is uniquely defined and 
continuous on [0, T] x R'^, the natural definition of a solution of the conservation law (1.2) with a given 
initial data mo belonging to the space of bounded measures A^i(R'^) is, see[16], 

m{t) = X{t;-)#mo. (2.2) 

Equation (2.2) means that the measure solution m{t) is the image (or the push-forward) of niQ by the 
flow X{t;-), i.e. for any Borel set B C M"', 

m{t){B) =mo{X{t--y\B)) =mo{{x eR'^ : X{t;x) £ B}) , 

or equivalently 

{m{t),^) = (mo, ^{X{t; •))) , for any (/) G Clj(R'') . 

It turns out that the definition of Filippov characteristics is well suited and that with it, (2.2) 
perfectly make sense as soon as the characteristics are unique. A definition of Filippov characteristics 
solution of (2.1) equivalent to the original one and easier to handle is given in [16], as follows: 
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Definition 2.2 (Filippov characteristics). For any given x and T > 0, a generalized solution 
the sense of Filippov of system (2.1) is a continuous map X{- ;x) : [0, T] i— ?• M'^ which satisfi' 

{M{a{;Vui-,t)) ■ v)}{X{t;x)) < dtX{t;x)-v< {M{a{-,Vu{-,t)) ■ v)}{X{t; x)) , a.e. t € [0, T] 

for any v ^W^, where 

\ M.{ai-, Vu(-.t)) ■ v)}(z) := sup I essinf a{y,Vu{y,t)) ■ v 

r>0 \y€B{z,r) 

{M{a{-,Vu{-,t)) ■ v)}{z) := inf ess sup a{y,'Vu{y,t)) -v]. 

''>0 \yeBiz,r) J 

Then, assuming 
(Ai) |a|GL-(M^;L-^p(M^)), 

the vector field a(-, Vu) satisfies (1.4) and the Filippov characteristics of system (2.1) exist. Moreover, 
if the following one-sided Lipschitz condition (OSL) holds true 

{a{x,Vu{x,t)) - a{y,Vu{y,t))) ■ {x - y) < -f{t)\x - a.e. t £ [0,T], x,y G M'^ (2.3) 

for a function 7 G L^{[0,T]), these characteristics are unique, the flow {t,x) 1— )■ X{t;x) is continuous 
on [0, T] X M'^, the map from R"' to M'^, x 1— )• X{t; x) is onto; the following existence and uniqueness 
result of a viscosity-measure solution of (1.1)-(1.2)-(1.3) can be stated: 

Theorem 2.3. Let uq G (VF^'~ n BUC){M.'^) and mo G Mi{'R'^). Assume (Hi)-(H2)-(H3) and 
(Ai). Then, if the viscosity solution u satisfies (2.3), there exists a unique measure solution m G 
C^{[0,T]; MiiR"^) - weak *) of (1.2)-(1.3) given by (2.2). 

In the reference case a{x,p) = p, as in (1.5), requiring the (OSL) condition (2.3) is equivalent to 
require that the viscosity solution u be semiconcave with respect to x, i.e. 

u{x + y,t)-2 u{x, t) + u{x -y,t) < (3{t) |yp , (2.4) 

for all x, y G M*^, t G [0, T] and for some /3 G L^([0, T]). The semiconcavity property (2.4) often 
characterizes the viscosity solution of Hamilton- Jacobi equations arising in control problems. However, 
for general vector fields a, it does not imply the necessary condition (2.3). Indeed, for Us = u * pe, 
where ps is the mollifier introduced above, the semiconcavity (2.4) implies: D'^u'^{x,t) v < P{t), for 
all V s.t. \u\ = 1 and for all {x^t) G Qt-, (see [7]). Moreover, Vu'^{x,t) —j' Vu{x,t) as e — ^ 0, for 
a. a. {x,t) G Qt- Owing to this convergence property, it is sufficient to obtain (2.3) for u^, whenever 
a is at least continuous w.r.t. p. Assuming in addition that a is differentiable w.r.t. p and satisfies a 
one sided Lipschitz condition w.r.t. x locally uniformly in p, we have 

a{x,p) - a{y, q) = Dpa{x, q + s{p- q)){p - q) ds + a{x, q) - a{y, q) , 
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and the estimate 



(a(x, Vn^(x, t)) - a{y, Vu'{y, t))) ■ {x - y) < 

ds da {x - yY V{s) D^u^y + a{x - y),t){x - y) + C \x - y\^ , 

with r(s) = Dpa{x,'Vu^ {y,t) + s{Vu'^ {x,t) — Vu'^ {y,t))). Therefore, the (OSL) condition follows easily 
from (2.4), if either Dpa{x,p) = I, as expected, or d = 1 and dpa{x,p) is non negative and upper 
bounded, as for (1.6). Apart from these two cases (already considered in [2]), the (OSL) condition is 
an assumption and has to be verified for the specific model at hand. 
We conclude this section by listing the properties that the measure solution (2.2) satisfies, see [16] : 

(i) monotonicity : mo > i^o =^ ^(t •) # rriQ > X{t ; •) ^ z/q; 

(ii) mass conservation : m(t)(M'^) = ?no(M'^); 

(iii) contraction property : \'m{t)\{B) < \mQ\{B) , for any Borel set B C W^] 

(iv) semi-group property : X{t - s ; •) # {X{s ; •) # mo) = X{t ; •) # mo; 

(v) uniform compactness at inGnity : V e > there exists R> s.t. |m(t)|(M'^\B/j(0)) < e, V i G [0, T]. 
It is worth to underline that the uniform compactness at infinity property (v) is fundamental to prove 
the convergence of approximate measure solutions toward the exact one in {[0,T]; A4i{M.'^) w — *). 



3 The semi-lagrangian scheme 

This section is devoted to the construction and the convergence analysis of a semi-lagrangian scheme 
for the time discretization of system (1.1)-(1.2) over the time interval [0,T]. For an introduction to 
this class of schemes we refer to [1, Appendix B] and [8, 9, 10, 11]. To proceed, we need to refine 
the previous hypothesis on the hamiltonian H and to add assumptions on the growth of H w.r.t. p. 
Therefore, let us assume in this section (Hi), (H3) and 

(H'2) for any R>0, supg^^^g^Q^^^ \H{x,t,p)\ = M{R) < +00; 

(H4) H is convex in p and either linear at infinity (i) or superlinear (ii): 

(i) 3 > s.t. lim ^^-f 'f''^^ = K, uniformly in {x^t) € Qt-, and there exists a positive 

|p|— i>oo 

function a € L°°{W^) such that \H{x,t,p) — H{y,t,p)\ < a{p)\x — y\ uniformly on x M'^ ; 

(ii) lim inf 'f'^^ = +00 . 

Remark 3.1 It is worth noticing that the hamiltonians (1.5) and (1.6) satisfy all the required assump- 
tions as soon as the potential V is uniformly continuous and bounded over Qt and Lipschitz continuous 
in X uniformly in t. The growth assumption (H4)-(i) can be relaxed to include also hamiltonians that 
are not uniformly positive at infinity. But we leave this generalization to the reader. 
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3.1 The semi-lagrangian scheme for the (HJ) equation 

Let E N be fixed and let h = T N^^ be the associated time step used in the semi-discrete scheme 
to be defined. Since the hamiltonian H is convex in p and continuous (in fact lower semi-continuity 
would be sufficient [17]), then H{x,t,p) = H**{x,t,p), where H* is the Legendre transform of H with 
respect top, i.e. H*{x,t,^) = suppgigdj^ • p — //(x, and the (HJ) equation can be written as 

dtu{x,t)= ini{-i-Vuix,t) + H*{x,t,i)}. (3.1) 

Next, plugging the first order forward finite difference for the approximation of dtu 

u{x,t + h) - u{x,t) 
dtu{x,t) 

and the first order approximation of the directional derivative ■ Vu 

^ „ , , uix — EhA) — uix,t) 
-i ■ Vu{x, t) ~ ^ ^ , (3.2) 

into (3.1), one easily obtains the following first order approximation for u{x,t + h) 

u(x,t + h) ~ inf {uix - ih,t) + hH*(x,t,^)} . 

Let us observe that whenever the exact solution u of (HJ) is semiconcave in x, then it possesses one- 
sided directional derivatives at any x and in any direction, i.e. the limit as /i \, of the right hand 
side of (3.2) always exists. Moreover, if u is differentiable w.r.t. x, the previous limit coincides with 
the left hand side of (3.2). In other word, the approximation (3.2) is consistent. 

For the semi-discrete scheme, it is sufficient to inductively define the approximation n" = u'^{x) of 
the exact solution u of (HJ) at time = n /i for n = 0, . . . , A^, by 

u''+\x)= inf {u''(x-^h) + hH*(x,e,n} , x € M^ (3.3) 

with = uq, the initial data for (HJ) in (1.3). Next, we need to show that it" as given by (3.3) 
shares the properties of uq, in the same way as the exact viscosity solution u does. These are well 
established facts in the context of the approximation of functions by the inf-convolution operator (see 
for example [7]). However, the difficulties here are on the one hand to show that the properties of 
the initial data uq are propagated to uniformly with respect to n and h, and on the other hand to 
handle the dependency of H on {x,t). 

Let us define Qh ■= x {0, . . . , A^} and the set of the arguments associated to the infimum in (3.3) 

A''{x):=arg mf {u''{x-Ch) + hH*{x,e,C)} , {x,n)£Qh. (3.4) 

Lemma 3.2 (Properties of H*). Let H satisfy (Hi), (H2) and (H3). // in addition H satisfies 
{H.4)-(i), then : 
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(a) H*{x,t,C) = +00 if \C\ > K, for any {x,t) G Qt , H*{-,-,0) G L°°{Qt) and H* is Lipschitz 
continuous in x uniformly in G B{Q,K) x [0, T]. 

On the other hand, if in addition H satisfies {il/i)-(ii), then : 

(h) H* also satisfies {H4)-(ii), for any r > there exists R = R{r) > such that 

H*{x,t,0= ma^ {p-^-H{x,t,p)}, V (x, t, G x 5(0, r) , (3.5) 

pGB{0,R) 

H* G L°°{Qt X i?(0, r)) and H* is Lipschitz continuous in x uniformly in (^, t) G -6(0, r) x [0, T]. 

The proof of the Lemma above is given in the Appendix. Let us just point out here that we have 
chosen to take the hamiltonian H to be Lipschitz continuous in x uniformly in p and t when H is hnear 
in p for \p\ — 7> +00, in order to obtain the Lipschitz continuity of H* w.r.t. x. Indeed, in this case the 
supremum defining H* is not a priori reached. Therefore we cannot conclude as in case (H4)-(ii). 

Lemma 3.3 (Properties of u"). Let uq G VF^'°°(M'^). Then, under hypothesis (Hi), (H2), (H3) and 
(H4), u" is well defined over by (3.3), i.e. A'^{x) is a non empty bounded set ofW^ uniformly in 
{x, n) G Qh and the infimum is a minimum. Moreover, G P^^'°°(M'^) and \\ u" ||^^i,oo(iud) is bounded 
uniformly in n and, assuming H* G L°°{Qt x B{0, K)) when the growth condition (H.4^)-{i) is satisfied, 
is Lipschitz continuous with respect to n uniformly in x, i.e. there exist three positive constants 
Co, Ci and C2 independent of h and n, s.t. 

||n"||^oo(Kd) < Co , ||Vn"||^oo(iRd) < Ci and (x) - {x)\ < C2\n - m\ h , (3.6) 

for all n,m = 0, ■ ■ ■ ,N, and x G M'^. 

Proof We suppose that tt" G W'^^°^{R'^) with 

II^"IIl°°(R'«) < C'o and ||Vn"||ioo(iRd) < C" , 

where Cq and Cf are independent of h. Then, the proof will follow by induction on n. 
It follows immediately by Lemma 3.2 that u^~^^ is upper bounded since 

'u"+^(x) < n"(x) + /iii'*(x,i",0) <CJ} + h\\H*{-,-,0)\\Lo.(^Q^), VxGM'^. 

Moreover, is obviously lower bounded since is bounded and H*{x,t,^) > —M, for any 

(x,t,^) G Qt X M'^, where M = supg^ \H{x,t,0)\, so that 

Ik"+^||L^(M.) <Co" + /imax{M, \\H*{;.,0)h^^Q^^} . 

Next, if H satisfies (H4)-(i), A"(x) C B{0,K) for any x G M'^ and the infimum in (3.3) is attained 
due to the continuity of u'^{x — ^/i) + /i i7*(x, t", ^) w.r.t. ^. 
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On the other hand, if if satisfies (H4)-(ii), since and H*(x, t, 0) are bounded and H* is superlinear, 
there exists > (increasing w.r.t. n but upper bounded uniformly in n and h) s.t. A'^(x) C 
-6(0, i?") for any x € M'^ and again the infimum in (3.3) is attained. 

Let us prove now that u""*"^ is Lipschitz continuous. We have for any x, y € and for a"{x) S A'^{x) 

(x) = (x - /i a" (x) ) + /i ii * (x, , q" (x) ) 

and 

n"+i(y) < u'^iy - ha'^ix)) + hH*{y,e,a''{x)) . 

Therefore, 

_ u^+\x) < [|lVn"||oo + hLH'] \x-y\< [C^ + /iL/f*] \x - y\ , 

where Lh* is the Lipschitz constant of H* w.r.t. x, and the statement follows exchanging the role of 
X and y. Consequently, € W^'°°{M.^). 

It remains to prove that is Lipschitz continuous with respect to n. As before, we have for any 
X G M'^ and the corresponding argument a"(x) G A^^x) 

|u"+i(x) - u"(x)| = [^"(x - /ia"(x)) - n"(x + /i i/*(x, t", a"(x))| 

< llV-u^lloo/ila^lx)! +/i|ii*(x,t",Q"(x))| < Cf/i|^"| + /i|ii*(x,t",a"(x))|, 

and the proof is completed, thanks to the uniform boundedness of A" and H* . □ 

Let us observe that when the hamiltonian H grows linearly at infinity, i.e. when the growth condition 
(H4)-(i) is satisfied, then H* is not necessarily upper bounded. Therefore, the additional hypothesis 
H* G L°°[Qt X B{{),K)) in Lemma 3.3 is necessary. It is easily seen that this additional hypothesis 
is satisfied by the hamiltonian (1.6). 

Finally, in order to obtain semiconcavity, we need the following semiconcavity hypothesis on H* . 



(H5) H* = H*{x,t,S,) is semiconcave in x uniformly in t and ^, with semiconcavity constant C, 



cone 



Lemma 3.4 (Semiconcavity). Let uq G W^'°°{M.'^) be semiconcave. Then, under hypothesis (Hi), 
(H2), (H3), (H4) and (H5), li" is also semiconcave uniformly in h and n. 

Proof. We proceed as in Lemma 3.3 supposing that u" is semiconcave with semiconcavity constant 
C^anc ■ Then, for any x G W^, any corresponding argument a"(x) G A^{x) and any y G W^, we have 
that 

n"+i(x + y) - 2n'^+i(x) + u''+\x - y) 

< vJ'ix + y- ha'^ix)) - 2u"(x - ha'^ix)) + u"(x - y - /la'^(x)) 

+ h [H*{x + y, t", a"(x)) - 2 H*{x, t", a"(x)) + ii*(x - y, t", a"(x))] 



< 



CZ-yrir -\- hC, 



H* 

cone 



\y\' 
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Hence, u""*" is semiconcave with C^}^]. < C^onc + ^^conc^ ^'^d by induction u" is semiconcave for all 
n = 1 , . . . , iV with : < + T C,^;, . 

□ 

Remark 3.5 Assumption (H5) is necessary since the hamiltonian H depends on x and t. It can be 
satisfied under regularity hypothesis on H. In the case of the hamiltonians (1.5) and (1.6), (H5) is 
satisfied if the potential V is convex in x uniformly in t. 

3.2 The semi-lagrangian scheme for the (CE) equation 

Now we can proceed to the construction of a semi-discrete scheme for the transport equation (1.2), 
replacing the time continuous fiow X{t;x) with a time discrete one. Set = vP' * and let us 
consider the following implicit backward Euler scheme 

=X" + /ia(X"+i,V<'+HX"+^)), n = 0,--- ,iV-l, 
= X , X . 

Next, given the initial measure mo and replacing (2.1) with (3.7), we define the semi-discrete ap- 
proximation of the measure solution (2.2) as 

m" = X"(-)#mo, forn = 0,...,iV, (3.8) 

i.e. for any Borel set B <ZW^ 

m'^iB) = mo((X")-i(S)) = mo{{x G M!^ : X"(x) G B}) , 

or equivalently 

(m",,^) = (mo,<^(X"(.))) for any € ^^(M'^). 

Above, the dependence of X" on x and e and consequently of m" on e has been skipped to simplify 
the notations, but it will be made explicit when necessary below. 

It is worth noticing that it is not possible to define the semi-discrete flow (3.7) with Vti""^"*^ instead 
of the gradient of the regularized u"^^, since we need X" to be defined for all (x, n) G Qh. Moreover, 
it is also not possible to use an explicit forward Euler scheme to define the trajectories X", because 
the discrete (OSL) condition (3.9) does not provide the necessary equicontinuity of these trajectories, 
in contrast with the implicit one (3.7), (see Lemma 3.6 and Remark 3.10 below). 

The following Lemma gives us the existence, uniqueness and the regularity properties of neces- 
sary to the well posedness of m" and to pass to the limit as h and e go to 0. Owing to the implicit 
nature of the Euler scheme (3.7), we have to assume an upper bound on the space-time step h. 

Lemma 3.6 (Properties of X"). Let uq G W'^'°°{W^) and assume (Hi), (H'2), (H3), (H4) and (Ai). 

Then, if in addition a G {C^{W^ x W^)Y , the (OSLh) condition 

{a{x, V<(x)) - a{y, V<(y))) ■ {x - y) < C\x - y\\ x, y G n = 0, • • • , iV , (3.9) 
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is satisfied, with constant C independent of e and h, and C h < 1, the solution of (3.7) is uni- 
vocally defined. Moreover, the flow {n,x) i-t- is locally bounded and Lipschitz continuous 

w.r.t. (x,n) G Qh, uniformly in e. 



Proof. The existence of X""^^ given X", is insured by the Brouwer fixed point theorem, the map 
y i-> X" + /ia(y, Vii"^^(y)) being continuous and bounded. The uniqueness of X^~^^ and the 
Lipschitz continuity w.r.t. x follow both from (3.9) and the upper bound on h. Indeed, for X"^^ and 
yn+i (Jgfing(i by (3.7) we have 

+ h (X"+^ - y^+i) • (a(X"+\ V<+i(X"+i)) - a(y"+\ v<+^(y'^+i))) 

< ix'^+i -y^+^iix" -y"| + c/i|x"+i (3.10) 

i.e. 

\ 1 — C h J 

Taking X" = we get the uniqueness, while iterating over n we get for two starting points x, y £ 

|X"(x) - X-{y)\ < (^1 + Y^Jlx - y\ . 

Therefore, for 6 G (0, 1) s.t. C h < 1 — 6, we have 

|X"(rE) -X"(?/)| <exp{nCh6~^)\x -y\ < eiip{CT 6~^)\x - y\ , n = 0,--- ,N. 

Finally, for any x G M"^ it holds true that 

\X"'{x) — X"^{x)\ < h\n — m\ sup sup \a{x,p)\, 

xeR'* pg"b{o,Ci) 

where Ci is given in (3.6), and we obtain the Lipschitz continuity w.r.t. n. The above estimate gives us 

also that the image by X^ of 5(0, R) is contained in the ball of radius {R + T sup sup |a(x,p)|). 

xeiR<* pe"B(o,Ci) 

Therefore, is locally bounded uniformly in e and n and the Lemma is proved. □ 

Lemma 3.7 (Properties of m"). Let mg G Aii{M.'^). Under the same hypothesis as in Lemma 3.6, 
the approximated solution in (3.8) is well defined in Moreover, m" satisfies the same 

properties (i)-(v) as the exact measure solution m, uniformly in n. 

Proof. By Lemma 3.6, the map x ^ X"'{x) is uniquely defined and continuous over W^. Moreover, 
it is onto from to W^. Therefore, the approximated solution m" is well defined in A4.\{W^). Fur- 
thermore, it is easy to see that m" satisfies all the properties (i)-(v) of the exact measure solution m. 
Concerning (v), it follows from the local and uniform in n boundedness of X" proved above (see 
Lemma 3.1 in [16]). □ 
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3.3 The convergence 

We can finally discuss the convergence of the semi-discrete scheme (3.3)-(3.7)-(3.8) to the solution of 
the continuous problem (1.1)-(1.2)-(1.3). The key tools to prove the convergence result have been given 
in [16], where the authors analyzed the stability of the measure solution with respect to perturbations 
of the vector field a. The stability of the Filippov characteristics is of course the basic stone. Later, 
these tools have been adapted in [2] to analyze the stability of the viscosity-measure solution of (1.1)- 
(1.2)-(1.3) with a{x,p) = p, with respect to the vanishing viscosity perturbation of the (HJ) equation. 
To discuss the convergence, we define the piecewise constant w.r.t. time approximated solution 

and Uy^ = Uh* Pe- Moreover, we also need to define the following time continuous trajectories by linear 
interpolation of and for any x G W^, 

X^(t;x) =X"(2;) + (t-r)a(X"+^(x),V<+^(X"+i(x))), te[t",t"+^], n = 0,---,A^. (3.11) 

Note that (3.11) gives us time continuous trajectories Xf^ with the same regularity as X", i.e. Xf^ is 
locally bounded and Lipschitz continuous w.r.t. (i,x) G [0, T] x M*^, uniformly in e and h. 

Theorem 3.8 (Convergence). Let uo G (Tyi'~ n BUC){W^) he semiconcave, mo G Mi{R'^) and 
assume (Hi), (H'g), (H3), (H4), (H5) and (Ai). Assume in addition that a G (C°(]R'^ x W^))'^, 
satisfies the (OSL^J condition (3.9) and Ch < 1. Let e = e{h) he s.t. e{h) — ?• as /i — )• 0. Then, 
as h ^ 0, 

(i) Ufi ^ u in L°°{Qt) and Vtt|(x,t) — )> Vn(x,t) at any point {x,t) G Qt of differentiahility of u, 
where u is the unique viscosity solution of (1.1); 

(a) Xf^ converges locally uniformly in Qt toward the unique Filippov characteristic X associated to 
the vector field a{-, Vu); 

(Hi) m in C^([0,T]; A4i(M.'^) w — *), where m is the unique measure solution of (1.2). 

Proof. Following standard results of the viscosity solution theory (see for instance [18]), it can be 
proved that there exists a constant C independent of h, s.t. 

-u(r)||ioc(]Rd) < n = 0,---,N. (3.12) 

Next, let {x, t) G Qt be a point of differentiability of u. For n = [|^] , by Taylor expansion and the 
semiconcavity of n" , we have that 

<(y) - <(x) - V<(x) ■{y-x)< {C^nc + ^^^c^Db - x?, V y G . 
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Therefore, since Vti"(x) p as h ^ 0, hy the uniform bound of Vu", the previous inequahty and 
the convergence of u^{t) toward u{t), gives us 

u{y, t) - u{x, t)-p-{y-x)< (C^nc + ^ C^oncW) \y - , 

i.e. p belongs to the subdifferential of u at x. Being u differentiable at {x,t), p = Vu{x,t) and we 
have obtained the proof of (i). 

We now prove the second part of the statement. From Lemma 3.6 it follows immediately that, up 
to a subsequence, Xf^ converges uniformly on every compact set of Qt to a continuous function Y, as 
/i — >■ 0. Next, since by (i) and the (OSL/j) condition (3.9), the (OSL) condition (2.3) is also satisfied, 
the Filippov characteristic X associated to the vector field a(-,Vn) is unique for any x € M'^. In 
order to prove that Y = X, let us define Af^{t;x) := Xf^{t;x) — X{t;x). We are going to prove that 
A|(f;x) — )■ as /i — >■ 0. Indeed, we have 

^dt\Al{t;x)\^ = dtXl{t-x) ■ Al{t;x)-dtX{t;x) ■ A|(i;x) . (3.13) 

By the Definition 2.2 of Filippov characteristics it holds true, for all x G M'^, a.e. t G [0, T] and all 
r > 0, that 

dtX{t; x) ■ Al{t- x) > ess inf a{y, Vu{y, t)) ■ A|(t; x) . 

yGB{X{t;x),r) 

Hence, for all 5 > 0, there exists x = x{5) € B{X{t;x),5) point of differentiability of u, s.t. 

dtX{t; x) ■ Al{t; x) > a{x, Vn(x, t)) ■ A^(t; x) - 5 . (3.14) 

Plugging (3.14) into (3.13) and using the definition (3.11) of Xf^, we obtain for = [^] ~ 1) 

i9i|A|(t; x)|2 < a(X"+i(x), V<+i(X"+i(x))) • A|(t; x) - a{x, Vu{x, t)) ■ A|(t; x) + 6 

= (a(X"+i(x),Vn^+i(^"+Ha;))) -a(x,V<+i(x))) • A|(t;x) (3.15) 
+ {a{x, V<+i(x)) - a{x, Vu(x, t))) ■ A|(i; x) + S := h + h + 6 . 

In order to estimate Ji, we decompose A|(t; x) as 

A|(t; x) = (Xf^it; x) - X"+i(x)) + (X"+i (x) - x) + (x - X{t; x)) , 
we note A := sup sup |a(x,p)|, where Ci is given in (3.6), and we make use of the (OSL^) 

xeB.'^ p£B{0,Ci) 

condition (3.9) to get 

h < 2A \Xf^it; x) - X"+i(x)| + C |X"+i(x) -x\'^ + 2A6<2A'^h + C |X"+i(x) - xp + 2A5 . (3.16) 

/,From (3.16), (3.15) and the uniform boundedness of A|(t; x) on every compact subset of Qt, we have 
obtained 

ldt\Al{t; x)\'^ <2A^h + C (x) -x\'^ + 2A6 + C|a(x, V<+i (x)) - a(x, Vu(x, t))\ + 5 . 

^1 
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Using the convergence results obtained above and passing to the hmit into the previous inequality as 
h ^ and 5 — > 0, we get that 



^dt\Y{t;x)-X{t;x)\ <C\Y{t;x) - X{t;x)\^ 



i.e. Y{t] •) = X{t] •) in L^^^(M'^) a.e. t G [0, T]. Finally, from the continuity of Y and X we deduce that 
Y = X. By the uniqueness of X and the uniform boundedness of Xf^, we have that all the sequence 
Xf^ converge. 

It remains to prove the convergence of m^. This can be obtained exactly as in Proposition 3.1 in [16] 
and we skip the proof. We just underline that the strong convergence of the trajectories Xf^, their 
time equicontinuity and the uniform compactness at infinity of m| (property (v)) are the key tools to 
prove that the sequence (m|(t), (f>) is convergent and equicontinuous on [0, T], for any € Cq(M'^). □ 

Proceeding as in (2.5), one can obtain the (OSL^) condition (3.9) for special a. This is summarized 
in the following Corollary which is a straightforward consequence of Theorem 3.8. It is worth noticing 
that the Hamiltonians (1.5) and (1.6) enter in the framework of this Corollary. 

Corollary 3.9. Assume the same hypothesis of Theorem 3.8, except the (OSLh) condition (3.9). If in 
addition, a is differentiahle w.r.t. p, satisfies a one sided Lipschitz condition w.r.t. x locally uniformly 
in p and either Dpa{x,p) = I or d = 1 and dpa{x,p) is nonnegative and upper bounded, then the same 
conclusions as in Theorem 3.8 hold true. 

Remark 3.10 (The explicit Euler scheme) The explicit Euler scheme 



does not provide a family of equicontinuous characteristics, under the (OSL/^) condition (3.9). This 
is a difficulty naturally intrinsic to the (OSL) condition that allows discontinuity of compressive type 
only, while the map x X'^{x) given by (3.17) can be expansive. Indeed, we have 



X' 



-n+l 



X" + /ia(X",Vn^(X")) 



(3.17) 



+ h (x"+i - y"+i) • (a(x", v<(x")) - o(y", v<(y"))) 



a(x", v<(x")) - a(y", v<(y"))p 



i.e., with the same constant A as in Theorem 3.8 



-n+l|2 



< (1 + 2C/i)|X" -y"|2 + 8>l2/j2 



and iterating over n : |X"(x) — X"(y)\ 



2 < g2CT [I 



x-y]"^ + 8A'^Th] . 
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4 The fully discrete semi-lagrangian scheme 



In this section we introduce a finite element discretization of (3.3) and an approximation of (3.8) by 
a bounded discrete measure, yielding a fully discrete scheme for (1.1)-(1.2)-(1.3). 

For an arbitrarily fixed space step k > 0, we consider the regular uniform grid of M'^ given by 
A"^ := {xi = ik , i G Z'^} . Let T'^ = {Sj}j^jk be the associated collection of non-degenerate, 
pairwise disjoint and uniform simplices having as vertices lattice points Xi G A**^ and covering W^, [Sj 
are triangles in dimension 2 and tetrahedra in dimension 3). We denote also by 



W'' = {we C{M.'^) : w is linear on S^, j G j''}, 
the space of continuous piecewise linear functions on T'^. Then, each w G can be expressed as 

Y,f^^{xMx^), (4.1) 



W[X) = 



for basis functions fi^ G satisfying Pfixj) = 6ij for i, j G Z'^. It immediately follows that any 
has compact support, < fi^ix) < 1, X^jg^d /3f (x) = 1 and at any x G M'^ at most {d + 1) functions 
13^ are non-zero. In the sequel, the interpolation operator defined by (4.1) will be denoted Pk- 

The fully discrete approximation of (3.3) based on the above space discretization is naturally given by 

f <(a^) = E.eZ'^/3f(^)<,, n = 0,--. ,iV, 

(4.2) 

ulf = ini^^^,{ul{x,-ih) + hH*{xi,t^,i)} , iGZ^ 

where ■ = uo{xi). It is obvious that, thanks to the continuous piecewise linear interpolation of the 
discrete approximation {u^^nzid on the lattice at any time step, the continuous function G 
shares the properties of the semi-discrete approximation vP' given in Lemma 3.3. Therefore, the 
equivalent of Lemma 3.3 for the fully discrete approximation will be skipped here. On the other 
hand, the semiconcavity of u^, given a semiconcave initial data uq, is not straightforward. Therefore, 
we shall give the equivalent of Lemma 3.4 here and the proof in the Appendix. 

Lemma 4.1 (Semiconcavity). Let uq G W^'°° (W^) he semiconcave. Then, under hypothesis (Hi), 
(H2), (H3), (H4) and (H5), vJj^ is discretely semiconcave for all n = 0, • • • N , i.e. 

ul{x + X,) -2ul{x) + ul{x - X,) < {C^Snc + TCZ,)\x,\^ , X G i G (4.3) 

and weakly semiconcave on W^, i.e. 



y\'^ + '^{E{x + y) + E{x-y)) 



ul{x + y)-2 ul{x) +ul{x-y)< (C7,"„V + T C^c) 
where E[x) is a nonnegative, continuous and bounded function vanishing overX^. 



x,2/G]R^ 



(4.4) 
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Remark 4.2 We have chosen to consider a uniform grid X , i.e. a uniform space step k in any axial 
direction, only for simplicity of notations. It is obvious that Lemma 4.1 still holds true if one chooses 
different space steps for each direction of the grid, but with an additional non-degeneracy condition. 
On the other hand, it seems difficult to obtain the discrete-semiconcavity of if the regular lattice 
is replaced with a general non-degenerate triangulation of M"'. Moreover, this key property strongly 
depends on the continuous piecewise linear interpolation (4.1). Therefore, it is not possible to use 
a nonlinear interpolation operator in order to preserve the semiconcavity of |xp and to obtain the 
semiconcavity of over from (4.3). 

We now turn to the approximation of (3.8). A definition of space continuous trajectories is however 
always necessary. Therefore, let ^ = * p;. and set 

= XI + ha{X^^\Vul^\X-^')) , n = 0, . . . , iV - 1, 

Again, the existence of XJ^~^^ given XJ^ is due to the Brouwer fixed point theorem applied to the map 
y I—)- XJ! + ha{y,\7u^~^^{y)). However, the same argument as in Lemma 3.6 giving the uniqueness 
of XJ^~^^ cannot be reproduced here. Indeed, due to the weak semiconcavity (4.4), ^ is only weak 
semiconcave in general. Therefore, we are allowed to assume only the following weak (OSL^) condition, 

{aix,Vul,ix))-a{y,Vul,{y))) ■{x-y)<C'\x-y\^ x,y GR^ \x - y\ > k , (4.6) 

for n = 0, • • • , and a constant C independent of h, k and e. 
In order to obtain the well posedeness of the implicit Euler scheme (4.5), we assume 

(A2) a locally Lipschitz w.r.t. p uniformly in x and one sided Lipschitz continuous w.r.t. x locally 
uniformly in p. 

Then, combining hypothesis (A2) with the Lipschitz property of n^, there exists a constant C", 
independent of h, k and e, such that if h is small enough, i.e. 

C"he~^<l, (4.7) 

the previous map is contracting, XJ^~^^ is unique and the flow (n,x) 1— )■ XJ^{x) is locally uniformly 
bounded and Lipschitz continuous w.r.t. (x, n) € Qh, uniformly in k. 

Remark 4.3 Condition (4.6) gives us nonetheless the interesting property that the characteristics X^ 
do not move much away from each other. Indeed, for Xj^~^^ and Yi^~^^ defined by (4.5) and such that 
- Y^~^^\ > k, proceeding as in (3.10), we obtain 

l^r^ - ^k^^\ ^ (1 + C'h6-^)\X]! - Y^\ . 

Therefore, 

l^n+i _ ^„,+i| < (1 ^ C'hd~')me.x{\X]: - y,"|, k} , 
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and iterating over n, 

\Xj:ix)-X^iy)\<exp{C'T6-^)max{\x-y\, k} . (4.8) 

Next, let 6i be the Dirac measure concentrated on the lattice points Xj. Let := Yli^i^f^^k i^i 
be an approximation of the initial measure mo in the space of discrete bounded measures, for the 
Mi{W^)w - * topology, conserving the mass and that tends to at infinity uniformly w.r.t. k 
sufficiently small, i.e. for all e > there exist > and K > d s.t. 

\ml\{R'^\BR{fd)) <e yk<K. (4.9) 

For example, one can consider ■ := mo{A^), with {A^)^^j^d a partition of such that Xj € iff 
i = j. Then, the mass is obviously conserved and (4.9) is satisfied since for any Borel set B 

\ml\iB) < \mo\ (Ui64^f ) > h = {i e Z'^ : x, € B} . 

We define /x^ as the image of by means of the flow XJ^{-). As for m", fi"^ is well defined in 
A^i(M'^), the map x i— > XJ^{x) being uniquely defined, continuous and onto from M.'^ to M"'. /x^ satisfies 
also all the properties (i)-(v) with respect to m^, for any n. However, //^ is not a discrete bounded 
measure. Hence, following the classical procedure of the finite element approximation, we observe first 
that to determine /x^, it is sufficient to test the equation 

{f,l<P) = {ml, <I>{X]: {.))), (4.10) 

against any (/> G C^(M'^). Indeed, C^(M'^) is dense in Cq (M'^) and /x^ tends to at infinity, uniformly in t 
and in k sufficiently large owing to (4.9). Moreover, since any function in C°(M'^) can be approximated 
uniformly by a function w € W'' with compact support, (4.10) becomes for such test functions 

= E E ml^^^iXJ:ixi))w{x,) = ml^w{x,) , (4.11) 

with rn^ ■ := '^i^Z'^^kif^ji-^ki^i))- observe that the last identity in (4.11) holds since the 

series X^iez-^ "^fc i absolutely convergent and for the same reason m^^- is well defined. Finally, (4.11) 
leads naturally to the following definition of the discrete bounded measure approximating fj,^ 

ml:=J2<iS^- (4-12) 
It immediately follows that conserves the mass, i.e. 

E <^ = E = "^o(^') ' n = 0, . . . , iV . 

Furthermore, property (4.9) is also satisfied by uniformly in n = 0, . . . , A^, by the definition of the 
coefficients m^- and (4.8). 
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With the above definitions, set UJ^ = {u^j)i(zid and = ("^^ JigZ'*- "^^^ fuhy discrete scheme 
(4.2), (4.12), reads 

C/,"+i= \ni{B\i)Ui: + hL^{i)}, n = 0,--- ,iV-l, 

(4.13) 

M," = A^M°, n = l, 

where i?'=(e) = (/3f (x, - K)^^^^,, L"(e) = (//*(xi, t", e)),6z<^ and = (/3f (Xnx,-))),jez- It is 
worth noticing that is a stochastic matrix, while A'^ is the transpose of a stochastic matrix. 
Defining again the fohowing time piecewise constant approximations 

Uh,k{x,t)=ul"'\x), {x,t)€QT, m^,(t) = mj^/^l , t e [0,T] , (4.14) 

as well as ^ = Uh,k * Pe and the time linear interpolated trajectories Xf^j^ exactly as in (3.11), we 
can prove the convergence of the fully discrete scheme (4.13). 

Theorem 4.4 (Convergence). Let uq £ {W^'°° n BUC){R'^) be semiconcave and mo G Mi{W^). 
Assume hypothesis (Hi), (H2), (H3), (H4), (H5), (Ai) and (A2). Assume in addition that the 
( 0SL\) condition (4.6) is satisfied with h sufficiently small such that C h < 1 and that (4.7) is 
satisfied with e = a G (0, 1). Then, as — ?> 0, 

(i) Uh,k u in L°°{Qt) and Vu^^(x,t) — ?> Vn(rc,t) in any point {x,t) € Qt of differentiability 
of u, where u is the unique viscosity solution of (1.1); 

(ii) Xf^ J. converges locally uniformly in Qt toward the unique Filippov characteristic X associated 
to the vector field a{-,Vu); 

(Hi) m'j^i^ m in C^{[0,T]; A4i{W^) w — *), where m is the unique measure solution of (1.2). 

Proof. We shall prove that there exists a constant C independent of h and /c, such that 

\\ul-u''\\L--{m.i)<C{n + l)k, n = 0,...,N. (4.15) 
As a consequence, by the time regularity of u and by estimate (3.12), it follows that 

\\uh,kii)-u(i)\\L^m <C{h + ^ + h^/^), vt e [o,r]. 

Estimate (4.15) is obvious for 7i = 0, being = PkUo and = uq. Let us suppose that it holds 
true for a given n. Then, for any argument a"(xj) € A'^{xi), we have by (3.3), (3.4) and (4.2) 

- 'i^^'^Hxi) < uUxi - h a"(xi)) - u'^ixi - h a"(xi)) < C{n + l)k . 

Exchanging the roles of u^^^ and u^~^^{xi) we obtain 

-u''^\xi)\<C{n + l)k, V^€Z^ 
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so that for any 2; € M'^, 

\ul+\x) - < |n"+i(x) - Pfcn"+i(x)| + ^ /3f - ulf\ <Ck + C{n + l)k . 

Since the convergence of Vu^ ^ can be proved exactly as in Theorem 3.8, from the weak-semiconcavity 
of ti^g, statement (i) is proved. 
Next, from (4.15), we have the same estimate for the regularized approximation u^^ and n", yielding 

\\VuI,-Vu^,\\l^^^,)<C^, n = 0,...,iV. (4.16) 

Then, from the (OSL^) condition (4.6), the (OSL/^) condition (3.9) with constant C follows as /c — )■ 0, 
and the trajectories X" can be defined as in (3.7), under the assumption C'h < 1. Moreover, 

+ /i(X^+i -X"+i) • (a(X^+\Vu^+i(X^+^)) -a(X"+\V<+i(X"+^)) 



+ /i(X^+i -X"+i) • (a(X^+\V<+^(X^+i)) -a(X"+\V<+H^"^^))) 
< - - X"| + /i(Lippa) - X"+1||Vm^+1(X^+^) - V<+i(X^+i)| 

i.e. using (4.16), 

|^n+l_^n+l| < (l + C'/j^-l) , 

for 6 G (0, 1) s.t. C' h < 1 — S. Iterating over n, we obtain for any x G M'^ 



|X^(a;) -X'^(x)| <C - ^{1 + C'h6~^y <C—^. (4.17) 

i=l 

Combining (4.17) with statement (ii) in Theorem 3.8, we have proved statement (ii). 

Finally, to obtain the convergence of ^ towards m, it is enough to prove that, as k ,h — )• 0, 
(m| /. — Trif^) ^ in C'^([0, T]; A4i{'K'^) w — *). The latter combined with statement (iii) in Theorem 3.8, 
gives us the claim. 

Again, using the compactness at infinity of (m| ^ — m^fj), uniformly in e ,h ,t and k sufficiently large, 
it is possible to consider only test functions (j) G C^(M"'). Then, for n = [|;] , we have 

{ml,it) - mlit), cP) = {ml - m", = {ml -m\4>- P,qy) + {ml - ?n", P^cP) 

= E - Pk<p{xi)) - (mo, 0(X"(.)) - Pkcl>{X'-{-))) + {ml - m", Pk<t>) 

:= h + h + h. 
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The term Ii is obviously equal to 0. The term I2 goes to as /c , /i — > by Lebesgue dominated 
convergence theorem. Concerning I^, using (4.11), we have 

/3 = {pi'^k^m - {m"',Pkct>) = {mlPu<P{Xl{-))) - (mo,Pfc0(X"(.))) 
= {mlP,cl^{X]i{-))-Pu<l^{X^{-))) + {ml-mo,Pk<t^{X^{-))) 
= E <^ miX^ix,)) - PkHX^'ix,))) + {ml - mo, Pfcc/>(X"(.))) := 4 + /^ 

By (4.17), the uniform continuity of Pkcf) and the uniform boundedness of m^, ^ as k ,h ^ 0. 
The same holds true for Jg since ^ niQ in A4i{'K'^)w — * as k ^ 0. To conclude, the sequence 
(m| ^(i) — mj^lt), (j)) converging to 0, is also equicontinuous on [0, T] since {mf^{t), 4>) is equicontinuous 
(see Theorem 3.8) and 

\{mlk{ti)-mlk{t2)A)\ < El"^M-<'N<^(^^)l- 

Using the uniform continuity of the basis functions, the time equicontinuity of the characteristics X^ 
and (4.17) in the estimate below 

\PKx''\x,))-P^{XI-{x,))\], 
we obtain the claim. □ 

5 Appendix 

Proof of Lemma 3.2. (Properties of H*). If if satisfies (H4)-(i), we need only to prove that //*(•, •, 0) G 
L°°{Qt), the other claims being straightforward. Then, let us define A = {p ^ M.'^ : H{x,t,p) < 
M , V (x, t) € Qt}, where M = supg^ \H{x, t,0)\. A is a non empty set of M*^ and it is bounded since 
H growths linearly for \p\ — > +00 and 

\P\ -H' 

Thus, there exists R> such that A C B{0,R) and 

H* {x, t, 0) = sup {-H{x,t,p)} = max {-H{x,t,p)} e L°°{Qt) ■ 

p£B{0,R) 



On the other end, if H satisfies (H4)-(ii), the uniform superlinearity of H* is a direct consequence. 
To prove (3.5), let us denote A{^) = {peR'^ : p-^- H{x,t,p) > -M , V {x,t) G Qt}. Then, A{^) 
is again non empty (0 G A(£^) for all ^) and bounded whenever |^| < r since 

H{x^ p-£M M 

— < 1 < r H , V p G j4 , 

\p\ \p\ \p\ \p\ 
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uniformly in {x,t) £ Qt- Finally, it follows immediately that H* € L°°(Qt x B{0,r)) and 

\H*{x,t,0-H*{y,t,0\<v{l + R)\x-y\, Vx,y GM^Vte [0,T],VeG5(0,r). 

□ 

Proof of Lemma 4-1- (Semiconcavity of u^). We follow here [8], [9] and [15]. The discrete semiconcav- 
ity (4.3) for x = Xj G X'' can be proved by induction. Indeed, this is true for by the semiconcavity 
of uq. Let a^- be one argument of the infimum in (4.2). Then, 

- 2 <j' + < uUx^+j - h al^) - 2ul{x, - h al,) + ul{x,., - h + h C,^;jx,f . (5.1) 

Next, since the point (xj — ha^-) belongs to a simplex Sj, it can be written as a convex combination 
of its vertices, i.e. (xj — ha'^j) = A^x;, where A; € [0, 1] and A; = 1. By the regularity of the 
lattice, the same holds true for {xi±j — ha^-) = Xixi±j. Finally, being piecewise linear on T^, 
(5.1) becomes 

I I I 

and (4.3) is proved on the lattice X^. For an arbitrary x € W^, we have again as before x = A/X/ 
and X ± Xj = Xixi±j. Hence, 

uUx + X,) - 2 n^(x) + ^^(x - X,) = E A/K+i- - 2 + , 
and (4.3) follows. 

Set C = {C^Snc~^f^Camc)^ inequality (4.3) implies that the continuous piecewise linear interpolation 
(^fc i ~ § l^«P)ieZ'* is concave on X^ for all n = 0, • • • iV, and consequently also over IR"'. Therefore, 
for any x, y G M'^ 

E pH^ + y)\^^\' - 2 E f^h^)\^^\' + E - 

Since the piecewise linear interpolation of [xp satisfies 

with £'(x) the approximation error, the weak-semiconcavity (4.4) is also proved. □ 
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ul+\x+y)-2ul+\x)+ul+\x-y) < ^ 
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